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Abstract 

Spiral galaxies are surrounded by a widely distributed hot coronal gas and seem to be fed by infalling 
clouds of neutral hydrogen gas with low metallicity and high velocities. We numerically study plasma waves 
produced by the collisions of these high- velocity clouds (HVCs) with the hot halo gas and with the gaseous 
disk. In particular, we tackle two problems numerically: 1) collisions of HVCs with the galactic halo gas 
and 2) the dispersion relations to obtain the phase and group velocities of plasma waves from the equations 
of plasma motion as well as further important physical characteristics such as magnetic tension force, gas 
pressure, etc. The obtained results allow us to understand the nature of MHD waves produced during the 
collisions in galactic media and lead to the suggestion that these waves can heat the ambient halo gas. These 
calculations are aiming at leading to a better understanding of dynamics and interaction of HVCs with the 
galactic halo and of the importance of MHD waves as a heating process of the halo gas. 
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1. Introduction 

Spiral galaxies are surrounded by hot coronal gas 
for whose existence three main origins are under de- 
bate: 1) Most plausibly, star-formation activity in 
the galactic disk (stellar wind and supernova type 
II explosions) pushes hot gas into the halo, by this, 
producing the hot-gas halo P, 0]. 2) Galaxy growth 
within the ACDM cosmology will inherently lead to 
gas accretion and by this produce a hot halo with 
virial temperatures if the cooling time is too long 
as for temperatures above 10^ K and thus for large 
galactic masses of M > lO^-'^AfQ. 3) And at least, 
cosmological simulations suggest that 30 - 40% of 
all baryons reside in a cosmic web of shock-heated 
warm-to-hot intergalactic medium accumulating as 
an extended hot halo around massive galaxies '3']. 

Another gas phase existing in the halo of spiral 
galaxies are high- velocity clouds (HVCs). In the 



Milky Way these gas cornplexes are detected in HI 
at high galactic latitudes [4| and are falling with high 
speeds towards the galactic disk. Since galactic halo 
gas is preferentially assumed to be produced through 
the above-mentioned process IV, some of this hot gas 
is suggested to condense, cool [5|, and to fall back to 
the galactic plane, or HVCs can stem directly from 
galactic HI gas that is swept-up by superbubbles, 
lifted into the halo, and returns on ballistic trajec- 
tories. While the first is, however, contradicted by 
Binney and collaborators because of heat conduc- 
tion, the second process would require also positive 
HI cloud velocities which are, however, not observed. 

This expected galactic fountain effect Q can, how- 
ever, not explain HVC velocities close or even above 
the galactic escape value. Moreover, HVCs contain 
clearly less than solar metal abundances and have 
probably larger distances [1, Q than achievable from 



Preprint submitted to Comp. Phys. Commun. 



January 13, 2013 



superbubble expansion, so that their origin from the 
galactic disk is very unlikely. In addition, there also 
exists a population of intermediate velocity HI clouds 
(IVCs). With their higher metallicity and far lower 
velocity of only 50 — 100 km ■ s~^ these clouds as 
well as those even slower (and smaller), so-called 
low- velocity clouds (LVCs), are also likely to reflect 
at least partly back-falling bullets expected from a 
galactic fountain. 

A not insignificant part of LVCs, however, can be 
assumed to have passed the galactic halo and is decel- 
erated by the drag of the ram pressure exerted by the 
hot halo gas [13] ■ This also leads to gas stripping from 
the clouds which is discernible as head-tail structure 
of HVCs As observed in several HVCs before 
and expected theoretically [H, the HVC 125-^41-207 
shows a two-phase structure with a low and high ve- 
locity dispersion for column densities larger than 2 x 
1015 m-2. 

Since in external galaxies seen edge-on like e.g. 
NGC 891 and NGC 4361 the hot halo gas is perme- 
ated with vertical magnetic field lines, and because 
the galactic hot halo gas is assumed to be mostly re- 
plenished by superbubbles which open perpendicular 
to the gaseous disk and, by this, also open the galac- 
tic magnetic field perpendicularly as it is coupled to 
the ionized gas, HVCs must be assumed to interact 
with the halo magnetic field. Unfortunately, topol- 
ogy, strength, and signatures of the galactic magnetic 
field are not well known and not yet manifested by 
observations. Consequently, the magnetic coupling of 
HVCs with the halo magnetic field is still vague and 
not yet elaborated. This innocence is not only caused 
by uncertainty of the field topology but mainly due 
to the total absence of magnetic field measurements 
in HVCs. 

Only Zimmer et al. jl3l | interpreted observations 
of X-ray gas connected with HVCs [14], and pref- 
erentially with the HVC complex C, as signatures 
of magnetic reconnection and explored the magnetic 
interaction of HVCs with the Reynolds layer 15|. 



Another and far more realistic interpretation of the 
multi-gas phase coincidence in HVCs is the sweep- 
up of hot halo gas by the large complexes on their 
passage th roug h the halo. HVC models by Vieser 
& Hensler in which the clouds survive the dis- 



ruption by Kelvin-Helmholtz instability due to self- 
gravitation, heat conduction, and cooling, the gas ac- 
cumulation from such surrounding hot gas has been 
demonstrated. 



Subsequently, Santillan et al. [17[ investigated the 
collision of HVCs with the magnetic galactic gas 
disk by means of 2D numerical models. In con- 
trast to the model assumptions of both studies no 
HVCs, i.e. large cloud complexes with velocities of 
—200 km • s^i and below, are found close to the galac- 
tic disk. Another limitation is that the magnetic 
field (in 2D) is only directed parallelly to the galactic 
plane. 

2. Motivation of the numerical study 

One of the unsolved problems in solar plasma 
physics is the heating of the solar corona. There exist 
mainly two possible explanations of this interesting 
problem - magnetic energy releasing and heating by 
the reconnection of magnetic field and heating the 
solar corona by plasma waves. In recent years, both 
of possible mechanisms were solved by many authors 
theoretically and also numerically (see e.g. 18 1, [l^ 
or fic' 




Figure 1: A HVC collides with the Reynolds layer, an ionized 
hydrogen layer of the galactic halo which includes magnetic 
fields. Picture taken and redrawn from |13|| . 



Since superbubbles expanding into the halo are 
supposed to cool adiabatically and the halo gas 
should continue at lower temperature to cool by radi- 
ation and by cloud evaporation, a continuous hot gas 
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supply is necessary. Because of the complexity of the 
multi-phase gas structure of the halo there is still a 
lack of sufficiently detailed quantitative explorations 
and results about its thermodynamics and magne- 
tohydrodynamics [21|. An additional heating source 
would therefore elegantly provide a solution to pre- 
vent the halo cooling. During the collisions of HVCs 
on trajectories oblique or perpendicular to the mag- 
netized galactic halo or galactic disk, the magnetic 
field lines along the galactic disc are distorted [13[, 



gas is compressed, and MHD waves are generated. 

Such MHD waves (including reconnection of mag- 
netic field) may, on one hand, lead to heating of the 
ambient gas and can, furthermore, influence the evo- 
lution and the dynamics of HVCs. This pro cess is 
analogous to the coronal heating in the sun [22| • This 
similarity of MHD processes encourage us to begin 
analogous MHD studies of HVCs passing through the 
halo and approaching the galactic disk while interact- 
ing with the magnetic field. 

3. Numerical model 

3.1. Simulations oj collisions of HVCs with galactic 
halo 

In plasma physics, there exist several methods, how 
the plasma dynamics can be described and calcu- 
lated, see e.g. ;235, {24]. 

For the calculations of HVCs' dynamics we used 
the full set of MHD equations [2^ with the addition 
of gravitational term: 



dg 
dt 



(1) 



dv 1 
e^ + g{v-V)v^~\/p+—{VxB)xB + gg, (2) 
dt ^0 



dB 

'dt 



V X (v X B), 



V5, 



dU 
'dt 
V -B = 



(3) 

(4) 
(5) 



Here p is a mass density, v (relative) flow velocity, 
p gas pressure, B is the magnetic held and g is the 



gravitational acceleration calculated by means of Eq. 
The plasma energy density U is given by: 



7-1 2 2/zo 



(6) 



with the adiabatic coefficient 7 = 5/3, and the flux 
vector S is expressed as: 

S=(u+p+^)-v^{vB) — . (7) 
V 2^0 / Mo 

The magnetohydrodynamic equations (HJ-Q were 
transformed into a flux conserving form, i.e.: 



a* dF{^) aG(*) 



dt dx 



dy 



0, 



(8) 



and were solved numerically. The vector ^ in our 
two-dimensional case is expressed as: 



/ P \ 

pVx 

pVy 

Bx 

By 

V u j 



(9) 



and G(^') are too com- 
for more information see 



The vector functions F(4') 
plex to be presented here, 
e.g. 0. 

For the numerical solution of this type of equations 
the two-step Lax-Wendroff algorithm is used. The 
numerical region is oriented in the x, y-plane, imple- 
mented at < a; < L and < y < L and covered by 
a uniform grid with 200 x 200 cells. We choose the x 
coordinate to be oriented in the galactic plane, while 
y points vertically into the halo. 

Open boundary conditions are applied and the 
time step satisfies Courant-Friedrichs-Levy condition 
in the form 1251: 



CFL • Aa 



max(cs -f- |v|) ' 
where the Courant number CFL is set to 0.8. 



(10) 
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3.2. Calculations of phase and group 

For the calculations of phase and group velocities 
of plasma waves in the galactic gas, we numerically 
solve the wave equation for plasma motions, where 
the equilibrium parameters (density and pressure) 
depend on y (e.g. j26i]): 



_d_ 

dy 



dy 



(11) 



Here Vy is the velocity component normal to the mag- 
netic field, w is the frequency, is the longitudi- 
nal wavenumber along the x-axis, and Alfven speed 

'^Aif = B^/i^og- 

The velocity Vx, parallel to the magnetic field is 
given by: 



dv„ 



(a;2 - klcl) dy 



(12) 



where Cg — {'^p/ qY^"^ is the sound speed. 

The function f{y) from Eq. ((TT]) is expressed as: 



/(2/) = 



Qcf{LO^ -klc\) 
(Cj2 _ klcl) 



(13) 



The tube speed ct and fast speed Cf are implied as 
CT = CsVMf/{cl + v\^if/'^ and Cf = (c^ + WAlf)^^^ 
respectively. 

Eq. pT|) has a singular point which is called the 
cusp resonance or cusp singularity. This point plays 
an important role in a case of slow magnetoacoustic 
waves, whereas it is not seen in numerical simulations 



of the fast magnetoacoustic waves (e.g. |27|). 

The second-order ordinary differential equation is 
rewritten in terms of two first-order equations of the 
new functions and 
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dy 



such that: 

de 



— = g{kxVj^if - uj 
dy 
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d6 
dy 



fivY 



The boundary conditions at the point y 
the "kink" mode are given by ^1 = 0,^2 



(14) 



(15) 

for 
c and 



the "sausage" mode satisfies the boundary conditions 
£,1 — c ■ /(0),^2 — 0, whereas the constant c is arbi- 
trary. 

To obtain a solution of Eq. ([Tl]) a fixed value of 
kx is used, integrating between y = and y = ?/max 
the two first-order differential equations and ([T5|) 
by means of the fourth-order Runge-Kutta method. 
Exact value of the frequency w is obtained by the 
bisection iteration method when the velocity Vy sat- 
isfies the boundary condition at the second point 
'^yiy = 2/max) = for both wave modes ("kink" and 
"sausage" mode). 

The total gas pressure ptot(gas and magnetic) is 
expressed as: 



-Ptot 



c?(a;^ — k^c^) dvy 



(w2 _ k^c'i) dy ' 



(16) 



and for the calculation of the magnetic tension force 
Tl we use the equation in the form, see (26j : 



lU! 



Q 



(17) 



This force is a part of the Lorentz force j x B and ap- 
pears whenever the magnetic field lines are curved. 
Magnetic tension force is determined by how much 
the magnetic pressure changes with distance (for de- 
tailed information see 26[). 



4. Initial conditions 

4-.1- Collisions of HVCs with galactic halo 

Numerical solutions of MHD equations are per- 
formed in 2D space. The absolute size of the simula- 
tion box amounts 3 kpc x 3 kpc. The initial position 
of modeled HVC is located 1250 pc above the galactic 
plane with an initial velocity of vo — —200 km • s~^ 

The mass density distribution of the galactic gas is 
given by the equation: 

g{y) = f?o[0.6e-«'/(280pc)=^o37g-yV(540pc)^_^ 
+ 0.1e-l^l/'*"°P'= + 0.03e-l^l/^"°P'=], (18) 
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and for the expression of the gravitational accelera- 
tion we use the equation [l7{ : 



9{y) 



8 • 10 
0.48e 



'^[1 -0.52e 

-|y|/900pci ^ 



-|a|/325pc 



(19) 



The galactic midplane gas density in Eq. p8|) is 
Qo — 2.24 X 10"^^ kg • m""^ and the initial magnetic 
field Bo = 5.0 /iG. The total pressure is given by the 
relation p{y) = J ggdy, with the boundary condition 
that Pout = p{y — 5.0 kpc) — Pa. The galactic mid- 
plane value of the total pressure po — 2.7 x 10^^^ Pa 

4-. 2. Phase and group speeds 

For both studied cases (uniform magnetic field par- 
allel to the galactic plane and galactic current sheet 
configuration) , we assume, in the state of equilibrium, 
initial plasma velocity v = 0. 

Equilibrium demands that the pressure (plasma 
plus magnetic) is uniform, i.e.: 



P 



El 
2^0 



= const. 



(20) 



4-.2.1. Uniform magnetic field 

In this case, the magnetic field is assumed B = 
Bq — const, in the whole simulation region. The total 
pressure p according to the Eq. (I20p is also assumed 
to be a constant, p ~ p^. 

4-2.2. Current sheet 

In the case of current sheet confi gura tion the mag- 
netic field is given by the equation [ISj: 



B = Bq tanh 



(y - V2) 



(21) 



and the equilibrium in Eq. (jl8p yields a plasma pres- 
sure given by: 



Piv) = Po sech^ 



(y - L/2) 



(22) 



given that p ^ (cold plasma) as {\y\/a) — > oo [26|, 
and where a = 0.5 kpc is the semi-width of the cur- 
rent sheet. 



Numerical solutions of Eq. PTI) are performed only 
for the so-called "kink" mode which well cor- 
responds to the situation of the HVC collision with 
galactic disk. 

5. Results and discussions 

In this section we present some selected numerical 
results obtained by means of our simulations. 

5.1. Collisions of HVCs with galactic halo 

In Figs. 2 and 3 we present the "time evolution" of 
a HVC collision with galactic media. Fig. 2 shows the 
perpendicular collision, whereas in Fig. 3 an oblique 
passage of the HVC with a = 45° through the galac- 
tic halo is depicted. 

The situations (a) show the gas density distribu- 
tion at the "beginning", i.e. Ta — 0.6 Myr af- 
ter the onset of the calculation, while the locations 
(b) present the density distributions for the time 
Tf, — 7.9 Myr in the case of perpendicular collision 
(Fig. 2) and Tb = 12.8 Myr for oblique cloud infall 
(Fig. 3). 



i 2 




Figure 2: The "time evolution" of the collision of HVCs with 
galactic halo in the case of perpendicular cloud infall. a: cloud 
at 0.6 Myr after the start of the models, b; position and struc- 
ture at time Ti, (see text). 



In both cases at time Tf, a dense regime in front of 
the HVC is clearly visible, where the galactic gas is 
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Figure 3: The "time evolution" of the collision of HVCs with 
galactic halo in the case of oblique cloud infall under a = 45° . 
a: cloud at 0.6 Myr after the start of the models, b: position 
and structure at time T;, (see text). 



1.5 3 
X [kpc] 

Figure 4: The total pressure in the numerical region and mag- 
netic field lines distorted by the HVC for the perpendicular 
infall of the cloud to the galactic disk at time Ti, (see text). 



compressed as the HVC moves towards denser galac- 
tic disk and is decelerated. 

Because this infall motion towards denser gas 
causes a bow shock to form, consequently, also 
plasma waves are created. 

In Figs. 4 and 5 are shown the total pressure of the 
galactic gas and magnetic field lines for both studied 
cases, i.e. perpendicular and oblique infall under an- 
gle a = 45°. The pressure on the head of the cloud is 
higher because the cloud moves through gas with in- 
creasing density near galactic disk, and the magnetic 
field lines are distorted because the magnetic field is 
"frozen" to the HVC. 

5.2. Phase and group speeds 

In Figs. 6 and 7 we present the phase and group 

speeds of wave signal, plasma velocity, total gas pres- 
sure and magnetic tension force for magnetic field ori- 
ented parallelly to the galactic plane and for a current 
sheet configuration, respectively. 

Notice, that in the graph of the phase and group 
speeds the scale on the left-hand side corresponds to 
the phase speed (solid line), whereas the scale on the 
right-hand side corresponds to the group speed (dash- 
dotted line) of the plasma waves. If we compare the 
phase and group speeds for both studied cases, we 




Figure 5: The total pressure in the numerical region and mag- 
netic field lines distorted by the HVC for the oblique infall 
(o = 45°) of the cloud to the gala<;tic disk at time Tj, (see 
text). 



can see that the profiles of phase speeds are similar, 
whereas the group speeds have different shapes. In 
the case of the current sheet the group speed becomes 
constant for high wavenumbers. 

From the figures of total pressure ptot we can dis- 
cern that for kxtt = 5 (red line) the lines are also 
the same, whereas for kxa = 10 (blue line) the total 
pressure is zero outside of current sheet, because of 




Figure 6: Phase (solid) and group (dash-dotted) velocity of 
wave signal (upper left) ; plasma velocity component Vy (upper 
right), total gas pressure (bottom left) and magnetic tension 
force (bottom right) for magnetic field oriented parallelly to 
the galactic plane, wavenumbers k^a = 5 (red) and k^a, = 10 
(blue). 




Figure 7: Phase (solid) and group (dash-dotted) velocity of 
wave signal (upper left); plasma velocity component Vy (up- 
per right), total gas pressure (bottom left) and magnetic ten- 
sion force (bottom right) for the current sheet configuration, 
wavenumbers k^a = 5 (red) and k^a = 10 (blue). 



constant (zero) plasma velocity Vy, see Eq. (jl6p . 

For the plasma velocity component Vy in Fig. 7 it 
is clearly visible, where the edge of the current sheet 



is located. Outside of the current sheet the plasma 
velocity is zero for kxO, = 10 and falls to zero in case of 
k^a = 5. This means that practically only the part 
of the galactic gas, where the magnetic field is not 
constant, is moving. In the case of a parallel magnetic 
field, galactic gas moves over the whole simulation 
region. It is probably caused by the fact, that inside 
the current sheet is present the Lorentz force j x B 
which does not allow the plasma waves generated in 
the center of the current sheet leave freely this space. 

As one can see, the y-component of the magnetic 
tension force Ti is in "anti-phase" to the plasma y- 
velocity component in all cases. 

6. Conclusions 

Our main aim is to investigate the effects of the 
collisions of HVCs with both the galactic media, hot 
halo gas and magnetic field, on their descendence to- 
wards the galactic disk. For this purpose here we 
present a first insight with respect to two aspects 
by 2D numerical models: at first, we solve numeri- 
cally the set of MHD equations and simulate the col- 
lisions of HVCs with the galactic media for two cases 
(oblique and perpendicular collision). In the second 
part of this paper, the equation of plasma motion is 
treated numerically in order to obtain the dispersion 
relations of waves in galactic media. 

In the literature [13| it is envisaged that the re- 
conncction of magnetic field lines in the galactic halo 
can lead to its gas heating. Unfortunately, but as a 
physically reasonable by-product, heating of the halo 
gas by plasma waves is not yet taken into account. 
There exists practically no paper dealing with this 
interesting problem. 

Our first results demonstrate that plasma waves 
provide an important heating process of the galac- 
tic halo gas in analogy to the solar corona. From 
the above-mentioned reasons our results are new and 
their further exploration of crucial importance for our 
understanding of the heating processes in the inter- 
stellar medium in general. Nevertheless, 2D simu- 
lations can only deal as a first insight but 3D sim- 
ulations are urgent for the MHD treatment to de- 
liver more accurate and reliable results. Therefore, 
as the next step we will advance our models to 3D 
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and to higher spatial resolution. For this purpose we 
intend to use as a powerful, well tested and widely 
applied tool for astrophysical simulations of many 
applications the FLASH code [2^, parallelized and 
with adaptive mesh refinement. Recently, [2^ per- 
formed first 3D models using FLASH, but with the 
main emphasis on the structural evolution of fountain 
clouds starting at a height of 5 kpc above the galac- 
tic disk. From the extension of our numerical models 
we expect important results and information about 
the MHD waves in galactic media. As further steps, 
also small-scale plasma processes like heat conduc- 
tion must be included, that leads to the stabilization 
of clouds fl6!| by suppression of KH instability. 
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